Measurement of the effect of quantum phase-slips in a Josephson junction chain 
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The interplay between superconductivity and 
Coulomb interactions has been studied for more 
than twenty years now 1 4 4 In low-dimensional 
systems, superconductivity degrades in the pres- 
ence of Coulomb repulsion: interactions tend to 
suppress fluctuations of charge, thereby increas- 
ing fluctuations of phase. This can lead to the 
occurrence of a superconducting-insulator transi- 
tion, as has been observed in thin superconduct- 
ing films ^, w ires^ and also in Josephson junc- 
tion arrays^EAHII]. The latter are very attractive 
systems as they enable a relatively easy control 
of the relevant energies involved in the competi- 
tion between superconductivity and Coulomb in- 
teractions. Josephson junction chains have been 
successfully used to create particular electromag- 
netic environments for the reduction of charge 
fluctuations^^. Recently, they have attracted 
interest as they could provide the basis for the re- 
alization of a new type of topologically protected 
qubit 17 18 or for the implementation of a new cur- 
rent standard 19 . Here we present measurements 
that show clearly the effect of quantum phase 
slips on the ground state of a Josephson junction 
chain. We tune in situ the strength of quantum 
phase fluctuations and obtain for the first time an 
excellent agreement with the tight-binding model 
initially proposed by Matveev et alP. 

The Hamiltonian for the theoretical description of su- 
perconducting circuits can be conveniently obtained by 
applying Devoret's circuit theory^. Here, each electrical 
element such as an inductance, a capacitor or the Joseph- 
son element can add a degree of freedom. In the case 
of circuits with a small number of electrical elements, a 
complete analytical description that takes into account 
all degrees of freedom can be obtained. However, when 
the circuits contain an increasing number of elements, as 
for example Josephson junction chains, even numerical 
solutions of the problem become difficult to obtain when 
taking into account all degrees of freedom. Nevertheless 
our measurements demonstrate that the ground state of a 
phase-biased Josephson junction chain (see Fig. [l|a)) can 
be described by a single degree of freedom. Although the 
chain is a multi-dimensional object, the effect of quan- 
tum phase-slips can be described by a single variable to, 
that counts the number of phase-slips in the chain. 

We start by giving a short introduction on the low- 
energy properties of a Josephson junction chain which 
have been studied in terms of quantum phase slips by 



Matveev et alP. Let us consider the Josephson junction 
chain depicted in Fig. [lja). The chain contains N junc- 
tions and is biased with a phase 7. We denote Ej the 

2 

Josephson energy of a single junction and Eq — fc its 
charging energy. Here we consider Ej ^ Ec- Let Qi be 
the charge on each junction and 0i the phase difference. 
In the nearest-neighbor-capacitance limit the Hamilto- 
nian can be written as: 

N N 

H = Y / l4E c (Q t /2e) 2 +E J (l-cos6 l )} ; ]T 0< - 7. (1) 

i=l i=l 

Ignoring the charging energy for the moment, we find 
the classical ground state, that satisfies the constraint 
on the phase Y^iLi^ii by minimizing the Josephson cou- 
pling energy. The corresponding phase configuration 
can be easily found for sufficiently large N (N 1) 
and it is given by Qi = j/N. The resulting Joseph- 
son energy hence reads E = -E/7 2 /2N and the chain 
is equivalent to a large inductance. If a phase-slip oc- 
curs on one of the junctions, say the jth junction, then 
9j — > 0j + 2ir. Since the Josephson energy is periodic in 
8j, a phase-slip of 2tt does not change the Josephson en- 
ergy of the junction j. However, the constraint Qi = 7 
is violated after such a phase-slip event. Therefore the 
phase difference Qi over all other junctions changes a lit- 
tle from j/N to (7 — 2ir)/N in order to accommodate 
the bias constraint. A phase-slip on a single junction 
leads to a collective response of all junctions. Conse- 
quently the Josephson energy of the entire chain changes 
from E = E, n 2 /2N to E 1 = Ej(>y - 2tt) 2 /2N. One 
can show similarly that the classical energy needed to 
accommodate to phase-slips without violating the con- 
straint is given by E m = Ej{^ — 2nm) 2 /2N . Therefore, 
the classical ground state energy of the chain, min {E m }, 
consists of shifted parabolas that correspond respectively 
to a fixed number m of phase-slips (see Fig. [ljb)) . For 
the special values 7 = 7t(2to + 1) the energies E m and 
E m+ i are degenerate. 

Taking now into account the finite charging energy Ec, 
quantum phase-slips can lift the degeneracy at the points 
7 = 7r(2m + 1). In the limit Ej ^ Ec, the hopping ele- 
ment for the quantum phase-slip can be approximated 
bjEHHl „ = l^E.jEclK{E,,/2E c ) 1/A e~^ E ->/ Ec . 
Since a phase-slip can take place on any of the N junc- 
tions, the hopping term between the two states |to) and 
|m + 1) is given by Nv. Therefore, using a tight-binding 
approximation, the total Hamiltonian for the chain is 



given by: 



ponentially suppressed with N and Ej/E(J^. 



H\m) = E m \m) - Nv [\m - 1) + \m + 1)] 



(2) 



Fig. [TJb) shows the numerical calculation of the two 
lowest eigenenergies of this Hamiltonian for three differ- 
ent ratios ^f- = 20, 3 and 1.3 in the case of a 6-junction 
chain. Fig. |t|c) shows the corresponding current-phase 
relation of the chain in the ground state. The chain's 
supercurrent is obtained by the calculation of the deriva- 
tive of the ground state energy E g : is — For 
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FIG. 1. (a) Schematic picture of the phase-biased Joseph- 
son junction chain, (b) Energy levels of a Josephson junction 
chain with iV = 6 as a function of bias phase 7 for differ- 
ent ratios Ej/Ec- For Ej/Ec = 20 (black lines) no split- 
ting is visible at the crossing points. For Ej/Ec = 3 (red 
lines) a gap emerges that increases rapidly with decreasing 
Ej/Ec- Blue lines show the energy levels for Ej/Ec = 1.3. 
For each Ej/Ec, the two lowest-lying states have been cal- 
culated by numerical diagonalization of the hamiltonian |2| . 
(c) Current-phase relation for the ground state £"3(7) for the 
same Ej/Ec ratios as in (a). The supercurrent is calculated 

from the derivative of the energy band: is = X "s -2 "' ^he 
chain current is reported in units of the critical current of a 
single chain junction io = t Ej. (d) Schematic picture of the 
chain shunted by the read-out junction, (e) Escape potential 
for the Josephson junction chain in parallel with the read-out 
junction for three different flux-biases <f>c in the read-out loop 
(see Fig. [2|. The ground-state of the chain clearly modifies 
the escape potential of the read-out junction. 

large values of Ej/Eq, quantum phase fluctuations are 
very small (v ~ 0) and the current-phase relation has a 
sawtooth-like dependence with a critical current that is 
approximatively N/tt times smaller than that of a single 
junction of the chain. We call this regime the "classical" 
phase-slip regime. When quantum phase fluctuations in- 
crease, i.e. Ej/Ec decreases, the current-phase relation 
becomes sinusoidal and the critical current becomes ex- 
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FIG. 2. Measurement circuit. The 6-SQUID chain is inserted 
in a superconducting loop. The flux $c created by on chip 
coils controls the phase difference 7 over the chain. The flux 
<&s through the SQUIDs can be controlled independently by a 
second coil. We denote the phase difference over the read-out 
junction S. 



Read-out junction 


SQUID at 4>s = 


S HU = (121±5)10 :i nm 2 


S bQ = (30 ± 2)lQ 6 nrri* 


C HU = 5.8T0.2/F 


C = 1.4T0.1/F 


Rf 3 = 968 ± 5fi 


R b N Q = 3800 ± 450^ 


I§° = 330 ± 2nA 


I% Q = 83 ± 9nA 



TABLE I. Parameters of the sample: size, capacitance, 
normal-state resistance and critical current of the read-out 
junction and a single SQUID of the chain. The critical cur- 
rent variance for the junctions in the chain is estimated to be 
smaller than 4%. 

To measure the effect of quantum phase-slips on the 
ground-state of a Josephson junction chain we have stud- 
ied a chain of 6 junctions. Our measurement setup and 
the junction parameters are presented in Fig. [2] and Ta- 
ble [T] Each junction in the chain is realized by a SQUID 
in order to enable tunable Josephson coupling Ej. In this 
way we can tune in situ the Ej/Eq ratio by applying a 
uniform magnetic flux $5 through all SQUIDs, and con- 
sequently we can control the strength of quantum phase 
fluctuations. For our measurements we placed this chain 
in a closed superconducting loop, threaded by the flux 
$c containing an additional shunt Josephson junction 
that is used for the read-out of the chain state. The flux 
$c enables the control of the bias phase 7 = $c — 5 over 
the chain. 

We have measured the switching current of the entire 
Josephson junction circuit containing both the chain and 
the read-out junction. The switching current was deter- 
mined from the switching probability at 50%. The escape 
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probability as a function of bias current hias shows the 
usual "S-form" with a width of w 20 nA. We apply typi- 
cally 10000 bias current pulses of amplitude Iu as and du- 
ration At < 1/is, and measure the switching probability 
as the ratio between the number of switching events and 
the total number of pulses. The results of the switching 
current measurements as a function of flux <S>c are shown 
in Fig. [3j From these switching current measurements we 
deduce the effect of quantum phase-slips on the ground 
state of the chain. 

The measured switching current corresponds to the es- 
cape process out of the total potential energy Utot con- 
taining the contributions of the read-out junction and the 
chain: 



of the measured switching current. Let us point out that 
we have used the nominal values for Ej and Ec calcu- 
lated from the characteristics of the sample indicated in 
Table |TJ The normal-state resistance for a single chain 
junction has been deduced from the measured normal- 
state resistance of the read-out junction by considering 
the size ratio between the two. We evaluate the precision 
of the determination of Ej and Ec to be in the range of 
±10%. This error bar on Ej and Eq yields an uncer- 
tainty SN = 1 for the effective number of junctions TV 
contributing to the phase-slip amplitude. This confirms 
the good homogeneity of our junctions and the collective 
nature of the phase-slip events. 



U tot (5,$ c ) =Ef°cos(S) 



E g (<S> C -d)- -iMasS. (3) 



Here E g is the ground state of the 6-SQUID chain calcu- 
lated by solving the Hamiltonian Q. As Ef° > E g the 
main component in Utot is the potential of the current- 
biased read-out junction Ej° cos(<5) — ^Ibiasb- Fig. ife) 
shows the escape potential at constant bias-currentfor 
three different flux values <pc corresponding to three dif- 
ferent biasing phases 7 over the chain. Let us point out 
that the position of the minimum of the potential U to t is 
in good approximation independent of the value of the 
flux <pc- Therefore the bias phase difference 7 over the 
chain depends only on the flux <pc ■ As a consequence, the 
4>c dependence of the measured switching current results 
from the 7 dependence of the chain's ground state. 

The escape from the potential U to t occurs via Macro- 
scopic Quantum Tunneling (MQT). The MQT rate 
for an arbitrary potential can be calculated in the 
limit of weak tunneling using the dilute instanton-gas 
approximation^. Within this model, the e scap e rate T 
out of the washboard potential U to t{l) reads^ 

T = Aexp[-B] , 
where A and B are given by: 



A=^™ 1 with 7 = 2/; 
c" r / gw 

JO Ly 16£g° U tot (x) 



B 



h 2 U tot {x) 

4 Eg° 



dx. 



dx 



(4) 



We have denoted by a the width of the barrier and by 
x the phase coordinate measured from the minimum of 
the washboard potential. The plasma frequency is ujq — 

\f 8E ° h tot ^ > where E^P is the charging energy of the 
read-out junction. 

Knowing the escape rate T, we can calculate the 
switching probability: 

P{hia S ) = i-exp [-r(4 



At]. 



(5) 



The results are shown as red lines in Fig. [3j The theory 
fits very well both in amplitude and shape the oscillations 
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FIG. 3. Measured switching current (black points) as a func- 
tion of (f>c over the chain for three different E.j/Ec ratios. 
The measurement noise for each point is about 0.2 nA. The 
red lines represent theoretical calculations for the switching 
current using formulas ([5| and Q. 

Fig. [4] shows the measured switching current ampli- 
tude AIsw and the corresponding theoretical calcula- 
tions as a function of Ej/Ec- For each measurement, 
Ej has been calculated using the flux dependence of the 
SQUID's Josephson coupling: Ej($ s ) = ^o{^s) with 

io(®s) = IcP cos J • In order to distinguish be- 
tween the suppression of the switching current that is 
due to quantum phase fluctuations and the one that is 
simply due to the well-known cancellation of the SQUID's 
critical current as a function of flux, we plot the switch- 
ing current amplitude divided by the critical current of 
a single SQUID iq. We see that the measured switching 
current amplitude follows very well the predicted theo- 
retical suppression of the switching current oscillations 
in the presence of quantum phase fluctuations. From our 
measurements we can also deduce the strength of the 
quantum phase-slip amplitude. With decreasing Ej/Eq 
ratio from 3 to 1 the quantum phase-slip amplitude in- 
creases from 0.8 GHz to 2.7 GHz. In addition, in Fig. 
[4] we have plotted for comparison the calculation for the 
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FIG. 4. Comparison between the measured (black points) 
and the calculated (red empty circles) switching current am- 
plitude as a function of the Ej/Ec ratio. Note that the 
switching current amplitude is divided by the flux dependent 
critical current of a single SQUID io, in order to reveal the 
effect of quantum-phase fluctuations. The top curve (blue 
empty circles) shows the theoretical calculation of the switch- 
ing current amplitude in the absence of quantum phase fluc- 
tuations. The lines are guides for the eye. 



switching current amplitude in the case when quantum 
phase fluctuations would be negligibly small: v ~ 0. As 
expected, we get a practically flat dependence on Ej/Eq- 

Further on, the upper x-axis of Fig. [4] shows the ra- 
tio Ej/ksT of the Josephson energy with respect to the 
thermal energy at T = 50mK. Since Ej ^> kgT, ther- 
mal fluctuations are excluded to explain the suppression 
of the switching current with decreasing Ej/Ec- Ad- 
ditional measurements (not shown here) reveal a con- 
stant switching current amplitude up to a temperature 
of T = lOOmK. 

In conclusion, we present for the first time a detailed 
experimental characterization of the effect of quantum 
phase-slips on the ground-state of a Josephson junction 
chain. These phase-slips are the result of fluctuations 
induced by the finite charging energy of each Josephson 
junction in the chain. The experimental results can be 
fitted in very good agreement by considering a simple 
tight-binding model for the phase-slips^. Our measure- 
ments also show that a JJ chain under phase bias con- 
straint can behave in a collective way very similar to a 
single macroscopic quantum object. 



These results open the way for the use of quantum 
phase-slips in JJ networks for the implementation 
of a new current standard, the observation of Bloch 
oscillations^, the fabrication of topologically protected 
qubits^ and the design of new superconducting circuit 
elements. 

Methods 

The circuit was fabricated on a Si/SiC>2 substrate and 
the Al/AlOz/Al junctions were obtained using standard 
shadow evaporation techniques. The aluminum oxide 
was obtained by natural oxidation in a controlled O2 at- 
mosphere. The sample was mounted in a closed copper 
block which was thermally connected to the cold plate of 
a dilution refrigerator at 50mK. All lines were strongly 
filtered by low-pass filters at the cryostat entrance and by 
thermocoaxial cables and tt filters at low temperatures. 

The switching current Isw of the circuit is obtained 
performing the following sequence. We use a series of M 
current steps of equal amplitude Iuas to bias the junc- 
tion. We count the number of transitions to the voltage 
state Msw and thus obtain the value of the switching 
probability Psw = A ^ w corresponding to the applied 
hias ■ By sweeping the Ibias amplitude and repeating the 
above sequence, we measure a complete switching his- 
togram, Psw vs hias- The Psw — 50% bias current is 
called the switching current of the circuit, Isw- 

The choice of the read-out junction critical current Iq° 
for an optimal measurement of is is not straightforward. 
On the one hand one would like to have Iq° ^> is, but 
on the other hand the width of the switching histograms 
w increases with 1^ and so does the probabilistic noise 
in the measurement w/\f~M. For reasonable measuring 
time scales, the number of current steps M is limited to 
values of about 10 . If we want to measure supercurrents 
for the SQUID chain in the range of InA, Iq° needs to be 
in the range of 100nA We have used a read-out junction 
with a critical current Iq = 330nA which offers a good 
trade off. 
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